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I.  INTRODUCTION 


Of  various  methods  proposed  for  removal  of  the  contaminant  droplets  on 
surfaces  of  a  vehicle  or  other  equipment,  utilization  of  water  jets  appears  to 
be  the  most  effective  and  the  most  economical  one  at  the  current  level  of 
technological  development.  The  procedure  is  to  use  the  water  jet  to  break  up 
the  droplet  and  subsequently  to  carry  away  the  contaminant.  The  surfaces  then 

•  can  be  decontaminated  by  moving  an  assemblage  of  water  jets  toward  the 

contaminant  droplets. 

_  In  order  to  design  an  efficient  jet  system,  the  following  knowledge  is 

vital:  the  general  flow  pattern,  the  evolution  of  the  contaminant  droplet, 

and  the  effect  of  each  flow  parameter  on  the  flow.  This  information  can  be 

sought  via  experiments,  but  that  would  require  very  sophisticated 

instrumentation  and  would  be  very  costly.  As  an  alternative,  one  may  use 
computer  simulations  based  on  appropriate  flow  models.  In  fact,  this  method 
can  have  a  much  greater  flexibility  than  experiments  for  examining  areas  of 
importance  in  the  flow  field  and,  therefore,  can  provide  better  insights  into 
the  flow  phenomena. 

The  jet-droplet  interaction  is  a  three-dimensional,  transient,  two-fluid 
flow.  The  flow  is  extremely  complicated  since  it  involves  interfaces 

separating  the  two  fluids  and  each  of  the  fluids  may  have  free  surfaces  with 
the  ambient.  In  the  present  investigation,  we  have  simplified  the  problem  by 
treating  it  as  a  two-dimensional  flow  for  which  we  have  developed  two  flow 
models,  namely,  one-fluid  flow  and  two-fluid  flow.  They  are  suitable  for  the 
computation  of  the  development  of  two  different  pre-impingement  flow 
configurations.  Results  presented  are  the  flow  pattern,  evolution  of  a 
contaminant  droplet,  effects  of  variation  of  flow  parameters  on  the  flow,  and 
some  typical  pressure  distributions  on  the  impingement  surface. 


II.  ROW  MODELS 


Figure  1  depicts  two  pre-impingement  flow  situations  which  can  occur  in 
the  decontamination  process.  In  the  first  configuration,  shown  in  Figure  la, 
a  water  jet  is  directed  at  a  contaminant  droplet  which  is  covered  by  a  water 
layer.  The  water  layer  is  stationary  or  flowing.  In  the  second  configuration, 
shown  in  Figure  lb,  there  is  no  water  coverage  on  the  contaminant  droplet.  To 
respectively  characterize  the  flows  developing  from  the  two  configurations  we 
have  developed  a  two-fluid  flow  model  and  a  one-fluid  flow  model.  Both  models 
describe  a  two-dimensional  viscous  flow  and  can  be  handled  by  the  computer 
code  SOLA-VOF  (Reference  1).  It  is  noted  that  the  computer  code  can  solve 
flow  problems  involving  fluids  separated  by  interfaces  in  a  region  without 
voids  or  one  fluid  with  voids  (ambient). 


1B.D.  Nicholas,  C.V.  Hirt,  and  R.S.  Hotchkiss,  "SOLA-VOF:  A  Solution 
Algorithm  for  Transient  Fluid  Flow  with  Multiple  Free  Boundaries,  "  Los  Alamos 
Scientific  Laboratory  Report  No.  LA-8355,  1980. 
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a.  Contaminant  Droplet  with  Initial 

Water  Layer  Coverage  (Two- Fluid  Flow) 


b.  Contaminant  Droplet  Without 
Initial  Water  Layer  Coverage 
(One-Fluid  Flow) 


Figure  1.  Pre-Impingement  Flow  Configuration 


A.  Two- Fluid  Flow  Model 


This  is  a  model  of  a  channel-type  flow  bounded  by  the  dashed  line 
indicated  in  Figure  la,  covering  the  major  part  of  the  flow  region.  The 
channel  contains  two  fluids  (water  and  contaminant)  separated  by  an 
Interface.  The  upper  wall  of  the  channel  coincides  with  the  upper  free 
surface  of  the  water  layer  so  as  to  eliminate  consideration  of  the  free 
surface  with  the  ambient.  An  outflow  boundary  condition  is  specified  at  this 
wall  and  at  the  ends  of  the  channel,  allowing  the  fluids  to  flow  out  the 
region*  The  contaminant  which  covers  a  rectangular  region  is  assumed  to  wet 
perfectly  the  lower  wall  of  the  channel.  To  account  for  viscous  effects,  a 
no-slip  condition  is  used  at  the  lower  wall.  Finally,  a  steady  uniform  jet 
velocity  is  specified  along  a  segment  of  the  upper  wall  as  shown  on  the  top  of 
Figure  2. 


B.  One-Fluid  Flow  Model 


Results  of  computations  with  the  two-fluid  flow  model  show  that  for  such 
a  close-in  impingement  the  interaction  flow  is  insensitive  to  the  variation  of 
viscosity  of  the  jet  fluid.  Also,  the  density  of  the  contaminant  (1070  kg/nr) 
is  very  close  to  that  of  plain  water  (1C00  kg/m3).  Therefore,  we  can  simplify 
the  problem  by  setting  the  physical  properties  of  the  jet  fluid  (water)  equal 
to  that  of  the  contaminant  and  treat  the  jet-droplet  flow  as  a  one-fluid 
flow.  The  flow  region  for  computation  is  bounded  by  the  dashed  line  in  Figure 
lb.  In  this  model,  only  the  free  surface  with  the  ambient  is  traced,  but  not 
the  Interface  between  the  water  and  the  droplet.  The  one-fluid  flow 

model  is  suitable  for  characterizing  the  flow  developing  from  the  configuration 
of  Figure  lb  for  which  the  two-fluid  flow  model  is  beyond  the  capability  of 
the  SOLA-VOF  code  because  the  model  requires  the  tracing  of  two  surfaces, 
l.e.,  the  interface  and  the  free  surface. 
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Figure  2,  Flow  Patterns  (Two- Fluid  Flow,  v  *  0.98  mm  /s, 
Vj  -  5  m/s.  e  -  45®)  W 

III.  BASIC  GOVERNING  EQUATIONS  AND  COHPUTATI ONAL  RESULTS 


The  basic  equations  governing  the  flow  are  the  two-dimensional,  unsteady, 
incompressible  Navier-Stokes  equations  and  the  continuity  equation.  They  are 
solved  in  the  SOLA-VOF  code  by  using  a  finite  difference  scheme.  The  free 
surface  (or  interface)  is  defined  by  a  function  F  which  satisfies  the  equation 
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where  u  and  v  are  the  velocity  components  in  x  and  y  directions,  respectively. 
This  equation  stai.es  that  F  moves  with  the  fluid.  The  value  of  F  is  unity  in 
computing  cells  fully  occupied  by  a  fluid  and  is  zero  in  cells  without  such  a 
fluid.  Cells  with  F  values  between  rero  and  one  contain  free  surfaces.  Using 
this  technique,  the  free  surface  can  be  tracked  by  monitoring  the  values  of 
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F.  In  addition,  the  code  allows  one  to  embed  Marker  Particles  in  the  region 
covered  by  the  contaminant  droplet,  providing  a  good  visualization  of  the 
evolution  of  the  droplet. 

The  dimensions  of  the  droplet  for  computation  were  assumed  to  be  3  mm  x 
0.6  ran  (length  by  height),  representing  the  average  size  of  the  droplet  on  a 
flat  surface.  Other  input  data  were: 


Dj  *  jet  size  *  1.83  mm 

V.  «  jet  velocity  “5-10  m/s,  corresponding  to  steady  dynamic 

pressure  of  14  -  50  kPa  (2  -  7.2  psi)  which  are  practical 
for  decontamination 


3  3 

p  “  contaminant  density  ■  1070  kg/m  ,  for  water  p  “  1000  kg/m 

c  w 

v  ■  kinematic  viscosity  of  contaminant  »  9.8  -  980 

mnr /s,  equivalent  to  10  -  1000  times  the  kinematic  viscosity 
of  water  v 

w 


IWo— Fluid  Flow  Model  (Contaminant  Droplet  with  Initial  Water  Layer 
Coverage ) 


Figure  2  shows  the  computational  domain  of  the  flow  field,  which  has  a 
dimension  of  20  mm  x  0.8  mm.  The  rectangle  in  the  flow  channel  is  the  region 
lnltally  covered  by  the  contaminant  droplet.  The  velocity  and  the  incidence 
angle  of  the  jet  are,  respectively,  V.  •  5  m/s  and  9  ■  45°.  It  is  noted  that 
all  of  the  graphs  in  this  figure  haveJ  been  magnified  by  a  factor  of  3  in  the 
vertical  direction  in  order  to  give  a  better  visualization  of  the  flow  near 
the  bottom  surface  of  the  channel  where  a  finer  computational  mesh  was  used. 
The  vectors  represent  the  local  velocities  in  the  individual  cells  of  the 
mesh.  The  figure  shows  the  flow  development  following  the  commencement  of  the 
jet  impingement  up  to  0.2  millisecond.  The  results  for  two  contaminant 
viscosities,  v  *  9.8  mmVs  and  u  ■  980  mm  Vs,  are  shown  in  the  left  and  the 
right  columns,  respectively.  Figure  3  provides  another  view  of  evolution  of 
the  droplet. 


The  viscosity  has  a  pronounced  effect  as  shown  by  the  droplet  profiles 
and  the  displacements,  S,  of  the  droplet  upstream  edge.  In  addition,  the 
fluid  particles  in  the  low  viscosity  droplet  have  been  disturbed  in  the  entire 
region,  while  in  the  high  viscosity  droplet  there  still  exists  a  region  near 
the  surface  in  which  the  fluid  particles  remain  in  good  order.  This  is  an 
indication  that  the  two  droplets  are  subjected  to  different  stress  levels. 


The  displacement  S  shown  in  Figure  3  is  defined  as  the  distance  from  the 
initial  location  of  the  upstream  edge  of  the  droplet  (marked  by  the  dashed 
line)  to  the  nearest  marker  particle  of  the  droplet  region  on  the  first  row 
above  the  impingement  surface.  We  embedded  24  rows  of  marker  particles 
uniformly  across  the  droplet  in  the  vertical  direction  and  there  were  30 
marker  particles  in  each  row.  As  illustrated  in  Figure  4,  the  nearest  marker 
particle  from  the  dashed  line  is  not  always  the  one  which  orginally  resided  at 
the  very  left  end  of  the  first  row.  This  is  because  the  jet  stream  is  lifted 
off  the  impingement  surface  and  the  marker  particles  in  the  front  are  carried 
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Figure  4 
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.  Movements  of  Marker  Particles  Fallowing  the  Start  of  the 
Jet  Flow  (Definition  of  the  Displacement  S) 
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upward.  The  location  of  each  marker  particle  following  the  jet  impingement 
can  be  traced  and  printed  out  by  the  computer.  The  displacement  S  defined,  in 
fact,  is  dependent  on  the  number  of  rows  of  marker  particles  embedded  and  the 
the  fineness  of  the  computational  mesh  size  used.  A  larger  number  of  rows  and 
a  finer  mesh  size  will  certainly  result  in  a  smaller  displacement  S.  In  the 
present  study,  due  to  the  excessive  costs  of  computer  time  and  the  problems  of 
numerical  instability  and  divergence  encountered,  the  finest  cell  size  in  the 
vertical  direction  used  near  the  impingement  surface  after  0.025  mm,  which  is 
approximately  4X  of  the  droplet  height.  The  initial  distance  between  the 
first  row  of  marker  particles  and  the  impingement  surface  was  also  0.025  mm. 
Test  runs  using  different  mesh  sizes  show  no  significant  changes  in  the 
droplet  profile  except  in  the  thin  layer  very  close  to  the  impingement 
surface. 

Figure  5  shows  the  effect  of  the  jet  incidence  angle  0  on  the 
displacement  S.  A  jet  at  a  smaller  incidence  angle  has  a  larger  portion  of 
the  jet  stream  moving  toward  the  droplet;  however,  it  produces  less  normal 

force  to  displace  the  droplet  because  of  a  smaller  rate  change  of  momentum  in 

the  impingement  arec.  A  jet  at  a  larger  incidence  angle  has  the  opposite 
result.  Therefore,  there  is  an  optimum  angle  at  which  a  jet  produces  the 
largest  displacement  S.  For  the  case  that  ■  9.8  mm^/s  and  Vj  *  5  m/s  and 
that  the  jet  is  fi.^ed  in  a  location  where  the  largest  S  can  be  obtained,  we 

have  found  that  a  jet  at  0  »  45  -  60°  gives  the  best  performance.  Further 

investigations  are  in  progress  for  cases  with  other  viscosities  and  jet 
velocities. 


B.  Pie-Fluid  Flow  Model  (Contaminant  Droplet  without  Initial  Water  Layer 
Coverage) 


We  next  consider  the  flow  configuration  of  Figure  lb  which  is  treated  by 
the  one-fluid  flow  model.  Figure  6  shows  a  sequence  of  flow  development 
following  the  initiation  of  the  impingement.  The  jet  stream  first  spreads  out 
on  the  wall,  then  engages  the  droplet  and  finally  is  lifted  off  the  wall  at 
some  angle.  The  figures  shows  free  surfaces  with  the  ambient,  but  no 
interfaces  between  the  jet  fluid  and  the  droplet  because  the  latter  interfaces 
are  not  traced  by  the  one-fluid  flow  model.  However,  the  evolution  of  the 
droplet  can  be  visualized  by  the  movement  of  the  marker  particles  shown  in 
Figure  7.  A  comparison  of  the  columns  in  Figure  7  shows  the  flow  dependence 
on  the  viscosity. 

Figure  8  is  a  plot  of  the  mean  velocity  S  of  the  droplet  upstream  edge 
versus  the  displacement  S  of  the  droplet  upstream  edge,  where  S  ■  S/t  and  t  is 
t'uc  time  for  the  displacement  reaching  to  S  after  the  start  of  the  jet  flow. 
The  Jet  velocity  strongly  affects  the  mean  velocity  S  of  the  droplet.  If  the 
viscosity  vc  "  vw  “  58  m^/s,  then  for  Vj  >  5  m/s  the  droplet  upstream  edge 
can  move  with  a  velocity  greater  than  3  m/s  at  its  early  movement,  aay,  in  the 
first  0.4  an  of  displacement.  The  velocity  is  large  enough  to  displace  the 
droplet  practically  instantaneously  at  the  impingement. 
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Figure  5.  Displacement  of  Oontamlnant  Droplet  Upstream 
Edge  vs.  Time  After  Commencement  of  Jet  Flow 


Figure  6.  Flow  Pattern  (One-Fluid  Flow,  v 
6  -  45°)  c 


9.8  mm  /b,  ■  10m 


Figure  7.  Evolution  of  Cbntaminant  Droplets 

(One-Fluid  Flow,  -  10  «/s,  6  -  45*) 


flL"A  *  «  *  J'j  irmr?  j? 


DISPLACEMENT  S,  mm 

Figure  8.  Mean  Velocity  of  Droplet  Upstream  Edge,  S,  vs.  Displacement 
of  Droplet  Upstream  Edge,  S 


Figures  9  and  10  present  a  comparison  of  flows  corresponding  to  two  jet 
sizes,  D,  *  D  and  D.  *  2.5  D,  for  the  jet  incidence  angle  0  »  0°.  We  see  that, 
the  additional  amount  of  the  jet  fluid  does  not  Increase  significantly  the 
displacement  of  the  lower  part  of  the  droplet  but  simply  flows  over  it. 
Investigations  for  cases  with  ©  >  0*  are  currently  underway.  Based  on  the 
above  results,  however,  we  can  expect  that  for  a  limited  jet  flow  rate, 
increasing  the  jet  velocity  rather  than  the  jet  size  will  give  a  better 
performance  In  decontamination. 

Another  area  of  concern  is  the  rise  of  the  pressure  peak  on  the 
lmplngaent  surface.  In  some  critical  areas,  such  as  the  optical  windows  of  a 
vehicle,  the  impact  pressure  that  the  areas  can  withstand  is  limited.  Figure 
11  presents  some  typical  pressure  distributions  on  the  impingement  surface 
after  the  start  of  the  jet  flow.  Because  the  phenomenon  is  transient,  the 
Instantaneous  pressure  peaks  can  rise  higher  than  jtwice  the  corresponding 
steady  dynamic  pressure  of  the  jet  velocity,  1/2  (pVj  ). 
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Figure  9.  Flow  Patterns  (One  Fluid  Flow, 

v  »  v  *  9.8  ma  /s,  V  *  10  m/s,  6  •*  0°) 
w  c  j 


Figure  10.  Evolution  of  Cqntaainant  Droplets  (One-Fluid  Flow, 
vc-  9.8  «r/s,  V  -  10  s/t,  0-0*) 
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Figure  11.  Instantaneous  Impact  Pressure  Distribution 
on  Impingement  Surface 

IV.  CONCLUSIONS 

The  computer  simulation  is  useful  for  observing  details  of  the  flow 
development  following  the  commencement  of  the  Jet  Impingement.  The  results 
indicate  that  the  contaminant  viscosit/  v£,  the  jet  velocity  V  ,  and  the  jet 
incidence  angle  8  have  strong  effects  on  the  profile  and  the  displacement  of 
the  contaminant  droplet.  The  computer  simulation  shows  that  for  V  -  5  m/s  and 
vc  ■  9.8  mm^/s,  a  Jet  at  an  incidence  angle  of  0  *  45  -  60“  gives  the  best 
performance  in  displacing  the  droplet.  Because  the  flow  is  transient  the 

Instantaneous  Impact  pressure  on  the  impingement  surface  may  rise  higher  than 
the  corresponding  steady-state  dynamic  pressure  of  the  jet.  Fbr  a  limited  jet 
flow  rate,  an  Increase  in  jet  velocity  rather  than  jet  size  will  likely  give  a 
better  decontamination  result. 

Farther  studies  on  surface  tension  effects,  moving  jets,  and  optimization 
of  flow  parameters  are  in  progress. 
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